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ABSTRACT 

This thesis provides a method to utilize existing digital 
simulations and missile reliability data to determine the 
probability of kill (Pk) versus range for the Seasparrow for each 
system state. A procedure is presented that uses the Pk versus 
range curve and the firing doctrine to determine the probability of 
defeating an antiship cruise missile attack for input into a 
baseline system effectiveness model. A method to optimize the 
firing doctrine using linear programming techniques is presented. 
Implemented using the General Algebraic Modeling System (GAMS), the 
optimized firing doctrine provides a significant enhancement to the 
probability of defeating an ASCM attack. The thesis enhances the 
Baseline Systems Effectiveness model for the NATO Seasparrow 
Missile System (NSSMS) developed by Professor W.M. Woods of the 


Naval Postgraduate School for the Naval Warfare Assessment Center, 
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THESIS DISCLAIMER 
The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that programs are free of computational and 
logic errors, they cannot be considered validated. Any application 
of these programs without additional verification is at the risk of 


the user. 
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I INTRODUCTION 


The NATO Seasparrow Missile System (NSSMS) is the primary 
point defense system for the United States Navy. Like many 
modern weapon systems, the Seasparrow is designed with a 
variety of alternate operating modes for use in the event of 
component failures. Although this is desirable from an 
operational standpoint, it creates difficulties in assessing 
the utility of the system. The Navy has traditionally tracked 
system availability by recording what percentage of the time 
a system is fully operational, usually referred to as A,. This 
method has several drawbacks when applied to the NSSMS. The 
major drawback is the lack of accounting for degraded modes of 
operations. Additionally, it provides no measure of the 
systems capability, only reliability and maintainability. To 
address these concerns, system effectiveness models were 
developed. In essence, a system effectiveness model sums, for 
all operating states, the product of the expected proportion 
of time a system will be in a given state (system state 
availability), the probability the system transitions to a 
different state during an engagement, and the capability of 


the system for the system state in which the system is 


operating at the time of missile launch. This paper provides 








an analytical framework to determine the probability of the 
NSSMS defeating an antiship cruise missile (ASCM) attack, for 
input into a NSSMS system effectiveness model developed for 
the Naval Warfare Assessment Center, Corona, CA (NWAC). This 
probability is one factor in the NSSMS system effectiveness 


equations. 


A. PURPOSE 

A method is developed to determine the probability of the 
NSSMS defeating an antiship cruise missile (ASCM) attack based 
on the operating state of the NSSMS at the time the missiles 
are launched. The thesis will explore the use of existing 
digitzl simulations developed by the Naval Air Warfare Center, 
China Lake, CA, (NAWC) to develop an estimate of the single 
shot probability of kill (Pk) for the Seasparrow when the 
NSSMS is in system state S; at the time of launch. This 
estimate will be combined with missile reliability data and 
used in conjunction with the firing doctrine to evaluate the 
overall probability of kill (Pk) for the NSSMS against a 
representative subsonic ASCM. Finally, a method to use linear 


programming techniques to optimize the firing doctrine versus 


a single target will be examined. 











II. BACKGROUND 

In an effort to improve the way in which the reliability 
of the NATO Seasparrow Missile System (NSSMS) is tracked, the 
Naval Warfare Assessment Center (NWAC) funded Professor W.M. 
Woods of the Naval Postgraduate School to develop a baseline 
systems effectiveness model for the system (Ref. 1]. The model 
he developed is based on methodology developed by the Weapon 
Systems Effectiveness Industrial Advisory Committee [Ref. 2]. 
The baseline model recognizes seven system states, S, to S,, 


and is written as 
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SE), »» AS;DS,, CS; (2.1) 
e1 jel 


where AS; is the operational availability of system state S;, 
DS;,; is the probability the system transitions from state S; 
to state S; during an engagement, and CS; is the capability of 
the system when in state S;. This paper provides a method to 
determine the value of CS; using digital simulation and 


missile reliability data. 


CS; is defined as the capability of the system to 
defeat the chosen threat given eight missiles are loaded at 
the beginning of the scenario [Ref. 1]. In the baseline model, 
CS; is defined as the probability of success given that two 


missiles are fired at the inbound target. This is written as 





CSy=P,(1) Pk, (1) +P, (2) Pk, (2) (2.2) 


where P;(n) is the probability of successfully firing n 
missiles and Pk;(n) is the probability of successfully 
destroying the target given that n missiles were fired in 
system state S,;. It should be noted that equation (2.2) 
incorporates the missile reliability term, R,, in the original 
cS; formula into the Pkj;(n) term for brevity. The formulation 
and calculation of CS; will be further discussed in Chapters 
IV and V. 
A. SYSTEM CONFIGURATION 

The NSSMS comes in two basic configurations: one fire 
control direc. or or two fire control directors. Each director 
has its own transmitter, receiver, signal data processor 
(SDP). and radar system controller (RSC). An eight cell box 
launcher is used in both configurations for launching the 
missiles. For purposes of systems effectiveness, if a two 
director configured ship has a casualty in only one director, 


the system is considered to be fully functional. 


B. SYSTEM STATES 

The baseline system effectiveness model for the NSSMS 
recognizes seven different system states. Although there are 
more than seven possible degradations to the system, only 


seven were considered to occur with sufficient frequency, or 


have a sufficient impact on probability of success, to justify 
inclusion in the model. The system states and their associated 
equipment casualties are listed in Table 1 [Ref. l,p. 6): 


Table 1: NSSMS OPERATING STATES 


System Operating Equipment Casualty 
State Mode 


All up except fault code 100: 
Transmitter power output degraded. 


All up except fault code 112: Range 
deviation too large. 


n 
un 


All up except FOC/RSC which has only 
co-.trol power sufficient to assign, 
tune, and fire missile in a normal 
manner . 


LI (Slaved All up except SDP or receiver. 
Tlluminator) 
or LLLTV 
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N 


All up except digital computer or 
complex Signal data converter (SDC). 
Casualty 
(ccc) 





Cc. MISSILE COMMANDS 

The RIM-7 Seasparrow missile is a derivative of the AIM-7 
Sparrow Air-to-Air missile. It has been optimized for the 
point defense role for use against antiship cruise missiles. 
The missile employs semi-active radar homing and pre-launch 
missile commands cannot be updated in flight. Degraded missile 
commands are the primary effect of equipment degradations in 


system states S,,S,,S,, and S,. The seven missile commands are 


listed below with one system command, launcher 
super~elevation, which is also affected by system degradations 


that may affect missile performance. 


a. Seeker Antenna (Head: Aim(pitch and yaw) 
b. English Bias 

c. Sweep Select 

d. Sweep Control 

e. Video Inject 

£. Range at Launch 

g. Channel Designate/HOH 


h. Launcher Super-Elevation 


Reference 3 provides a complete discussion of these missile 
commands. In the event of a system casualty that prevents the 
fire control system from providing inputs to the missile, 
state S, for example, default values are utilized. Of these 
commands, Range at Launch is evaluated as having the greatest 


effect on missile performance. 


D. DIGITAL SIMULATIONS 

Due to the limited number of actual firings, it is 
necessary to conduct simulations to determine the missiles 
probability of success. The NATO Seasparrow program office 
(NSPO), maintains a network of simulations which are available 


throughout the Seasparrow technical community [Ref. 4]. The 


Navy organization responsible for the Seasparrow simulations 
is the Naval Air Warfare Center, China Lake, CA (NAWC). NAWC 
maintains two Seasparrow flight simulations that offer 
sufficiently accurate modelling to give valid results: a six 
degree of freedom digital simulation and a six degree of 
freedom, hardware-in-the-loop simulation. Due to the order of 
magnitude different in cost between the digital simulation and 
the hardware-in-the-loop simulation, the digital simulation is 
the preferred choice. In addition to the flight simulation, 
NAWC maintains several endgame lethality simulations that will 
be utilized, including a target detecting device (TDD) 
simulation, a warhead fragmentation pattern simulation, anda 
target vulnerability model. 

The four simulations are used sequentially to determine 
the probability of success for each simulation trial. The 
flight simulation is used to determine the flight path of the 
missile from launch to intercept. The last few seconds of the 
missile and target flight paths are then entered into the TDD 
simulation to calculate the warhead detonation point. After 
determining the warhead detonation point, the relative 
positions and velocities of the missile and target are fed 
into the warhead simulation to determine how many fragments 
struck the target. The estimated fragmentation damage to the 


target is then compared with the target vulnerability model to 


determine the probability of the target being destroyed. 





IIt. SIMULATION PLAN 

Due to the limited number of actual missile firings versus 
ASCMs, simulations are required to determine the impact of 
range and system degradations on the probability of kill (Pk) 
for the Seasparrow in each system state. Digital simulations 
offer the advantage of substantially lower cost and 
Significantly shorter turn around time, when compared to 
hardware-in-the-loop simulations. 

In any systems effectiveness model, a specific scenario is 
required to determine the probability of success of the 
system. The baseline systems effectiveness model [Ref. 1) 
provides three Seasparrow missiles configurations and three 
possible ASCMs to be evaluated. The simulation plan in this 
thesis is developed for the RIM-7P configuration of the 
Seasparrow versus the Exocet cruise missile. This selection 
was chosen since all elements required (flight simulation, TDD 
simulation, warhead simulation, and target vulnerability 


model) are readily available at NAWC China Lake [Ref. 5). 


A. BACKGROUND 

Initially, it was assumed that the probability of kill for 
the Seasparrow could be represented by a single value for the 
entire range of the missile. Conversations with the technical 


community [Refs. 5&6) indicated that this would be an over 


simplification. This would be especially true of any degraded 


modes of operation for the system. It is expected that Pk will 


increase as the velocity of the Seasparrow increases, and 
decrease when the velocity of the missile decreases. This 
effect should be especially noticeable at or near maximum 
range when the missile is coasting and has very little energy 
available for terminal maneuvers. It is anticipated that the 


Pk versus range curve will look roughly as shown in Figure 1. 


Figure 1: Expected PK vs Range Curve 


It should be noted the values for Pk indicated above are 
for illustrative purposes only and do not represent the actual 


probability of kill for the Seasparrow missile. 





B. REQUIREMENTS 

The simulation plan must conform to several constraints. 
The plan must cover the full range of the Seasparrow, have a 
sufficient number of test points to provide good shape 
definition for the PK versus Range curve, have enough trials 
at each point to provide a reasonable confidence interval for 


Pk, and meet the cost constraints discussed in Reference 1. 


C. ASSUMPTIONS 

As a result of limitations of the digital simulations, 
several assumptions were made when the simulation plan was 
created. First, the engagement is in an electronic counter- 
measures (ECM) clear environment, at least in the spectrum 
covered by the NSSMS director. The engagement is limited to a 
single inbound ASCM, and each missile is assumed to act 
independently from all others. Finally, the firing ship is the 
target for the ASCM, and the Exocet flies its optimum profile, 
detecting the firing ship at maximum range and commences sea 


skimming as soon as possible. 


D. SYSTEM STATES 

Due to the scenario chosen, RIM-7P versus Exocet in a ECM 
clear environment, the seven system states can be reduced to 
five states by combining states S,, S, and S,;. Degradations to 
the FOC, S,, are only important in an ECM environment, while 


a computer or signal data converter failure, S,, only affects 


10 


the ability to acquire the target and should have negligible 


effect on the probability of success. 


E. SIMULATION PLAN 

To explore the Pk versus range dependency, thirteen 
intercept points were chosen, with the points spaced closest 
at the extremes of range. Intercepts closer than five percent 
of the maximum range are not evaluated as this is inside the 
range of the close in weapons systems. A sample size of fifty 
was adopted as a compromise between cost and confidence limit 
precision. Complete details of the Simulation Plan can be 


found in Appendix A. 





IV. DATA ANALYSIS 


Once the simulations are completed, the resultant data and 
the missile's reliability data needs to be analyzed and 
combined to determine the Pk versus range values. Once these 
values are determined, they can be compared to the firing 
doctrine to find Pk(n),;, the probability of killing the target 
given n missiles are successfully fired in system state Sj, 
and P,(n), the probability of successfully firing n missiles 


during an engagement. 


A. ESTIMATING PROBABILITY OF KILL 


An estimate of the probability of kill for any single 
missile is given by combining the probability that the missile 
guides correctly once it launches, Ry, the probability that 
the missile doesn't fail in flight, Ry, and the probability 
of success given that first two conditions are met, Pk,. The 
reliability information can be derived from missile flight 
data available from the flight analysis department at NWAC, 
while the probability of success is derived from the 
simulations. 

1. Determining Pk, 

The fifty values at each simulation point constitute 


fifty independently, identically distributed (IID) estimates, 





bkg;, i=1,2,...,50, of the probability of kill, Pk,. By the 
Central Limit theorem, the mean, x, of the fifty estimates 
should have a nearly normal distribution with mean Pk,. If o 
denotes the standard deviation of the probability distribution 


for Pks, then the quantity 


ime Pe): (4.1) 
(a//n) 

has a nearly standard normal distribution. Since the standard 

deviation, o, is not known, the sample standard deviation, s, 

is substituted for go in equation (4.1). The probability 

distribution of the resulting expression is approximated by a 


t-distribution with n-1 degrees of freedom. The resulting 


formula for a 95% lower confidence limit, ks,,), for Pk, is 


Bk,,,, =x-—t,, ; 4.2 
'S(L) J ra 0.95 ( ) 
where 
— 1 2 
ar S Pk,;, (4.3) 
ce Eas 
rece > (Bk,,-X)?, (4.4) 


and n is the number of trials (replications of simulations) 
and tin.1),0.95 is the 95th percentile of a t probability 


distribution with n-1 degrees of freedom . 
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2. Missile Reliability 

After Pk, is determined, it must be multiplied by the 
two missile reliability terms to determine the value for Pk at 
any intercept range. The first term, the probability that the 
missile guides once it is launched, Ry,, is simply the 
probability that all components of the missile function 
properly at launch. The second term, the probability of no 
in-flight failures, Ry, varies with the duration of the 
flight [Ref. 7]. Missile reliability can be adequately 
represented by an exponential life distribution reliability 
model, where the probability of a successful flight is defined 


as 


P(Reliable flight) =Ry Rye (4.5) 


where Ry, is the probability that all missile components are 
functioning correctly at launch and Ry is the probability 
that a missile functions correctly until intercept. Ry has an 
exponential life distribution. 

a. Estimating Ry, 

A point estimate, Ry, of Ryo the probability of the 
missile to correctly guide once launched, is the ratio of 
guidance successes divided by the total number of launches. A 
guidance failure is defined as a missile that launched but 
failed to guide properly because not all components functioned 


at launch. The number of successes is the number of tests 


14 





minus the number of guidance failures. A 100(1-a@)% lower 
confidence limit, Rygi,), for Ryg can be obtained by looking up 
the appropriate value in a table of confidence limits for 
binomial attribute data using existing tables such as 
Reference 8. Alternatively, 


n-f 


ee -.) (4.6) 
(n-f) + (£+1) Fy ge (2 (n-£+1) ,af) 


Regs tz) = 


where n is the number of launches, f is the number of 
failures, and F,., is the 100(1-a)% percentile point of the 
F-distribution with 2(n-f+1) and 2f degrees of freedom. If 
RIM-7M and RIM-7P missiles have sufficient commonality of 
components, the number of test firings could be pooled for 
determining Ryg and Rygiz)- 
b. Estimating Ry 
The second component of missile reliability, Ry, 
involves failures that occur in flight. Due to the limited 
number of failures that could be attributed to in-flight 
failures, as opposed to undiscovered failures that occurred 
prior to launch, a qualitative statement on the distribution 
of in-flight failures is impossible. For simplicity, the 
failure rate, 4, is assumed to be constant, which implies an 
exponential life distribution for the in-flight time to 


failure. The resultant term for in-flight reliability is 


Rygze* (4.7) 


and its estimate Ry is 


5 








a (4.3) 


t is the time of flight to the intercept point, f;, is the 

number of in-flight failures, and T; is the total flight time 

for all missile test firings. Again, if the RIM-7M and RIM-7P 

have sufficient commonality, the test results can be pooled. 
3. Estimating Pk, 

The estimate of the probability of kill derived from 
the simulations is multiplied by the estimate of the 
probability that the missile guides correctly, and the 
estimate of the missiles in-flight reliability for each of the 


thirteen sample points; that is 


Bk, =RyRePky, 1=1,2,...,13 (4.10) 
where Ry is recalculated for each intercept point. The result 
is an estimate for Pk at each of the thirteen intercept 
ranges. It is then necessary to fit a curve to the points to 
allow estimation of Pk at any possible intercept range. 

4. Curve Fitting 
There are variety of commercially available graphics 
programs that perform curve fitting. Given the expected shape 


of the Pk versus range curve, a high order, seventh to tenth, 
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polynomial curve should best fit the data. Figure 2 below 





shows a ninth order polynomial curve fit generated by the 







cor mercially available program DeltaGraph Professional, but 







many other graphics programs support curve fitting. 






Pk vs Intercept Range 
Ninth Order Polynomial Curve Fit 
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Figure 2: Sample Curve Fitting 








When the best fit has been determined for the data, 
the formula used to generate the curve provides an estimate 


for Pk for any intercept range. 


B. ESTIMATING Ry, 
Once the value of bk have been calculated for all possible 


intercept ranges, the next component required to calculate cS; 





is P,(n), the probability of being able to launch n missiles. 
The factor required to calculate P,(n) is Ry, the probability 
that the missile successfully launches. A success is defined 
as a missile that, after the firing signal was generated, 
launched. A failure is when the launch signal is generated, 
and the missile does not leave the launch rail. Since the 
NSSMS detects misfires and automatically attempts to fire the 
next available missile, these failures should not 
significantly affect the probability of kill, but do impact 
the probability of being able to launch a specified number « : 
missiles. 

A point estimate, Ry, of Ry, the probability of the 
missile to correctly launch, is the ratio of the number of 
successful launches divided by the total number of attempted 
launches. A 100(1-a@)% lower confidence limit, Ryip), for Ry, 
can be obtained by using tables, as for Rygi,)- Alternatively, 


ee ec (4.11) 
Ng, + (Map-Mgp +1) Fy ig (air (61) 2 iaqp-ge) ) 


where ng, is the number of successful launches, ng, is the 
number of attempted launches, and F,_, is the 100(1~a)% 
percentile point of the F-distribution with 2(n,,+1) and 
2(nay-Ng,) degrees of freedom. As with the other reliability 
terms, if the RIM-~7M and RIM-7P missiles have sufficient 
commonality of components, the number of attempted firings 


could be pooled for determining Ry, and Ry (1) - 
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Vv. CAPABILITY OF SYSTEM 


Once the estimates for the value for Pk at any intercept 
range, Pk;, and the reliability of the missile to launch, Ry, 
have been made, it is possible to estimate CS;, the capability 
of the system in state S;. In the baseline model CS; was 
defined as the probability of success given two missiles were 
fired. Since the NSSMS has sufficient range to allow a shoot- 
shoot-look-shoot-shoot doctrine against most threats, CS; has 
been redefined as the probability of success given up to four 
missiles are fired during an engagement. This change allows 
CS; to more accurately reflect the expected employment of the 
system versus subsonic cruise missiles. CS; is made up of two 
components, Pk(n),;, the probability of success given that n 
missiles are fired in state S;, and P,(n), the probability of 
successfully firing n missiles given eight missiles are 
available at the beginning of the engagement. The new equation 


for cs; is written as 


4 
CS5=>> P, (2) Pk(n) 4, (5.1) 


nel 


where n is the number of missiles fired during the engagement. 





A. ESTIMATING P, (n) 

To estimate the probability of firing one, two, three, or 
four missiles, the missile launch reliability term, Ry, is 
required. For the case of one, two, or three missiles fired, 
the binomial formula provides the probability of each 


occurrence. 
B(n) = 7 )Re (2- )9", ns3. (5.2) 


Pi, (4) includes the probability that four or more missiles are 


able to be fired, and is written as 


£,(4)-5° (8)Ag(1-8_)**. 


uw 


To simplify the calculations, 6,(4) can be thought of as one 
minus the probability that n is less than four, which is 


written as 


3 
B, (4) =1-( ($B, (n)) + (1-Re) ®) « (5.4) 
nel 


Once each term for PL(n) has been estimated, it is multiplied 


by the appropriate Pk(n), term to find Cs;. 








B. ESTIMATING Pk(n),; 

To estimate the value of Pk(n);, the firing doctrine is 
reviewed to determine at what range the intercepts would 
occur, and the appropriate values of hk; are determined from 


the Pk versus range curve. The formula for &k(n) 5 is then 


Pk(n) =1-T] (1- (Bk,|RI,) ) (5.5) 
ei 


with RI, being the estimated range at intercept for missile i. 
Once the value of fk(n); has been determined for all values of 
n, the values can be multiplied by their respective f,(n) 
terms and summed to find CS;. When estimating k(n); each 
missile is assumed to act independently, and their is no 
accumulation of damage to the target. This simplifies the 
calculation and provides the most conservative estimate of 
success. 

For system state S,, SI/LLLTV mode, on those ships that do 
not have a three dimensional radar to provide positioning data 
to the NSSMS director, replenishment ships for example, the 
low light level television is used for tracking, and the 
maximum range of detection, and thus intercept, must be 
shortened accordingly. Specifications for LLLTV expected 


detection ranges can be found in Reference 9. 
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C. SEARCH RADARS 

The detection range plays an important role in determining 
the capability of the system. This is especially true for the 
Computer Complex Casualty mode, state S,, which relies on 
manual acquisition of the target with the NSSMS fire control 
radar. Adding search radar performance and acquisition times 
to the baseline system effectiveness model would enhance the 
accuracy of the baseline model, but would entail a substantial 
increase in the complexity of the model. Although not covered 
in detail, a method to incorporate search radar performance is 
provided as a basis for future enhancements to the baseline 
system effectiveness model. 

1. Modifications to Baseline Model 

The baseline system effectiveness model is modified to 

account for search radar performance and reliability by 
introducing an additional factor in each term of the original 


formula. The system effectiveness equation becomes 


m #7 7 
SE). p> > AC,AS DS, 3CS5, x (5.7) 
1 g°1 Jel 


where AC, is the availability of each shipboard search radar 
configuration C,, m is the number of possible configurations 
of the search radars, and CS; , is the capability to defeat 


the threat given that the NSSMS is in state S;, and the search 





radars are in configuration C,. The other terms remain 
unchanged from the baseline model. 

Each search radar configuration, C,, is determined by 
the search radar configuration of the ship. A Spruance class 
destroyer might have eight configurations as shown in table 2. 
A Wichita class replenishment oiler, AOR, might have only two 
configurations, MK23 Target Acquisition System (TAS) up and 


TAS down. 


Table 2: DD963 SEARCH RADAR CONFIGURATIONS 


Radar Status 


TAS up, SPS-40 up, SPQ-9 up 
TAS down, SPS-40 up, SPQ-9 up 


TAS up, SPS-40 down, SPQ-9 up 


TAS down, SPS~40 down, SPQ-9 up 


TAS up, SPS-40 up, SPQ-3 down 
TAS down, SPS-40 up, SPQ-93 down 


TAS up, SPS-40 down, SPQ-9 down 





TAS down, SPS-40 down, SPQ-9 down 


Depending on the search radar in question, there could 
be several possible operating states that could effect the 
capability of the NSSMS to engage the target. For instance, on 
some ships, the TAS computer links the combat direction system 
and the NSSMS together [Refs. 9&10). If the TAS computer 


fails, the NSSMS operators must manually acquire all targets, 
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increasing the acquisition time significantly. This would 
result in three states for TAS, all up, radar down but 
computer up, and all down, instead of two, and the number of 
configurations would increase accordingly. 


The new term for the capability of the system, CS; ,, 


needs to incorporate three additional components; the expected 
range of first detection, the acquisition time, and the 
expected probability of success given the target is acquired 
at a specific range. The addition of these terms would allow 
the effects of system states S,, FOC/RSC degraded, and S;, CCC 
mode, to be isolated from system state S,, all up, resulting 


in a better ectimate of the capability of the system. 
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VI. OPTIMIZING FIRING DOCTRINE 


Once the Pk versus range curve is known, it is possible to 
optimize the firing doctrine to maximize the probability of 
defeating the threat. Ideally, if the optimization were a 
linear program, one could take advantage of commercially 
available optimization software such as GAMS or LINDO. Pk(n) 
is maximized by minimizing the probability of missing. The 


probability of missing can be written as 
4 
IT (1-(Pk,|RI,)) . (6.1) 
#1 


Although the probability of missing is still a non-linear 
equation, by taking its natural logarithm, it becomes a linear 


equation, 
4 
De tn (a (Pky|RT,)) , (6.2) 
=i 


which is readily solvable as a modified set covering problem. 


A. LINEAR PROGRAM 
The first requirement to allow the use of linear 
programming techniques is to break up the Pk vs range curve 


into discrete time intervals with constant Pk within the 


interval. The use of one second intervals based on the time of 





first possible launch (maximum range intercept) until the last 
possible launch (minimum range intercept) results in 
appreximacely one hundred and twenty time intervals. This 
provides sufficient resolution to minimize the effects of 
discretizing the continuous data, but makes the problem small 
enough to be readily solvable. For each launch time interval, 
the appropriate intercept range as the percentage of maximum 
range is entered into the formula selected for the curve 
fitting to return the Pk for the time interval. Solving the 
following formulation results in the optimum firing doctrine. 
1. Index 
t: The time interval in which the missile is fired. t is 
indexed from 1, the maximum range intercept launch, to 
tmax» the minimum range intercept launch. 


u: An alias for t, u offers the same index for use with 
matrix notation. 


2- Data 


P.: The natural log of the probability of missing (1-Pk,) 
for a missile fired during time interval t. 


Cy,.: A confounding matrix, a one indicates a missile 
launched at time interval u will be airborne during t= 
interval t. A zero indicate a missile launched at time 
interval u is not airborne during time period t. An 
additional five seconds is added to each time of flight 
to allow for operator assessment of the engagement. 


3. Decision Variable 


X,: A binary variable, a one indicates a missile should be 
launched during time interval t, a zero indicates no 
missile is to be fired during time interval t. Although 
X, is a binary variable, in practice it only needs to be 
limited to values less than or equal to one to provide 
binary results. 
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4. Formulation 


Caan 
Minimize p> PX, (Obj) 
t= 


Subject to: 


Coax 
Dy * <4 (1) 
c= 
od 
> Cu,%, s 2, for all u (2) 
tel 
X,-. + X, < 1, 2stst,,, (3) 
X,s 1, for all t (4) 


In these equations, constraint (1) limits the number of 
missiles fired to a maximum of 4 per engagement, constraint 
(2) limits the number of missiles airborne at any time to two, 
constraint (3) limits the rate of fire to one missile every 
two seconds, and constraint (4) limits the number of missiles 


fired during any given time period to one. 


3B. RESULTS 
To determine the effect of using the optimized firing 
doctrine, the linear program was implemented using the General 


Algebraic Modeling System (GAMS). To avoid using classified 





data, a range of ten nautical miles (NM) and a constant 
velocity of 600 nautical miles per hour (kts) was assumed for 
the Seasparrow. The Pk versus Range curve found in figures 1 
and 2 was utilized, and the Exocet was assumed to have a 
constant speed of 600 kts. The optimized firing doctrine 
returned a value for Pk(4); of 0.967. A shoot-shoot-look- 
shoot-shoot doctrine with the first shots at near maximum 
range returns a value of 0.927, a difference of 0.04. An even 
bigger difference is seen when the probability of destroying 
the target with the first two missiles is examined. With a 
shoot when the target is within range doctrine, the 
probability of destroying the target is 0.584. The optimized 
doctrine increases the probability of success to 0.752. This 
is especially significant given the limited number of missiles 
available in the launcher and the extended time required to 
reload the launcher. In all cases, the increased Pk from the 
optimized doctrine results from a shorter intercept range, 
with a corresponding increase in single round Pk. 

Although these numbers are only for illustrative purposes, 
they show the significant increase available from an optimized 
firing doctrine when compared to a doctrine that assumes a 
constant Pk. The optimization program ran on a 80486/33MHZ 
computer in just under 1.5 seconds, and can be easily modified 
to accommodate almost any firing doctrine by simply changing 
constants. The GAMS model is provided in Appendix B. Before 


using with actual data, the missile in-flight confounding 
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matrix, C(u,t), needs to be updated with actual time of flight 


data. 





VII. CONCLUSIONS AND RECOMMENDATIONS 


Dr. Woods' baseline system effectiveness model for the 
NSSMS provides the framework to assess both the availability 
and capability of the system. This thesis provides an 
enhancement to the baseline model by providing a method to 
determine the probability of kill for the NATO Seasparrow 
missile as a function of range and system state. The Pk versus 
range curve can be developed using existing missile 
simulations and, when used in conjunction with the available 
missile reliability data, can provide a way to estimate the 
capability of the NSSMS to defeat the designated threat. 

Beyond enhancing the system effectiveness model, a Pk 
versus range curve can provide insight into the proper 
tactical employment of the system. Through the application of 
linear programming techniques, it is possible to significantly 
increase the probability of success by changing the firing 
doctrine to maximize Pk. Additionally, the development of Pk 
versus range curves for degraded states of operation can 
provide the basis for firing strategies to minimize the 
effects of the degradations. 

Further enhancements to the baseline system effectiveness 
model should concentrate on incorporating the effects of 


search radars and acquisition times into the model, as 
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outlined in chapter V. The baseline model could easily be 
expanded to include NSSMS, search radars, and close in weapons 
systems (CIWS), providing an assessment of the capability of 
a ship to defend itself against a chosen threat. 

Modern weapon systems with several degraded modes of 
operation make the traditional methods of tracking systems 
availability less suitable. The baseline systems effectiveness 
model for the NSSMS, with the enhancements provided by this 
paper, provides the technical community with a better method 


of measuring the utility of the NATO Seasparrow. 
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APPENDIX A. SIMULATION PLAN 
The following simulation plan is proposed for execution by the 
Naval Air Warfare Center, China Lake, California, using the six 
degree of freedom, digital Sparrow simulation and associated 
endgame lethality simulations. The procedure for requesting 
utilization of simulations in the NSPO performance assessment 


models network to conduct studies is outlined in Reference 1. 


Common to all cases. 

Missile: RIM-7P 

Target: Exocet SSM 

Target Profile: Exocet profile with firing ship as the intended 
target. The Exocet acquires the firing ship at maximum range 


and commences sea skimming at the earliest opportunity. 


Simulation Points: 


5.0 50 

7.5 50 
10.0 50 
15.0 50 
20.0 50 
30.0 50 
40.0 50 
55.0 50 
70.0 50 
85.0 50 
90.0 50 
95.0 50 
100.0 50 


All range at launch missile command degradations are 
created by adjusting the actual launch range, not the range at 
launch commands. The range at launch command should always reflect 
the range required for intercept at the desired range. 
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Case 


Case 


Case 


Case 


Case 





1: NSSMS system states S,: All up, S,: FOC/RSC Degraded, and 
S;: Computer Complex Casualty Mode. 


Missile Commands: Nominal 
Launcher Super-Elevation: In accordance with algorithm 
Transmitter Power: Nominal 


2: NSSMS system state S,: Transmitter power degraded. 


Missile Commands: Nominal 
Launcher Super-Elevation: In accordance with algorithm 
Transmitter Power: Degraded 


3: NSSMS system state S;: SDP range signal degraded. 


Missile Commands: Range at launch deviation +6%-15% from 
actual, all others nominal 

Launcher Super-Elevation: In accordance with algorithm 

Transmitter Power: Nominal 


4: NSSMS system state S,: SDP Range deviation too large. 


Missile Commands: Range at launch deviation greater than 
+15% from actual, all others nominal 

Launcher Super-Elevation: In accordance with algorithm 

Transmitter Power: Nominal 


5: NSSMS system state S,: Slaved Illuminator/Low Light 
Level Television mode. 
Missile Commands: Default 


Launcher Super-Elevation: Default 
Transmitter Power: Nominal 


Note: Values for missile command degradations presented here 


and in Reference 1, Appendix F are tentative, and will require 
modification following the NATO Seasparrow technical community 
review of Reference 1. 





33 


APPENDIX B. GAMS SOURCE CODE FOR OPTIMIZED DOCTRINE 


STITLE Thesis Formulation 
SSTITLE Optimized Firing Doctrine 
*------------- -GAMS AND DOLLAR CONTROL OPTIONS --------------------- 


SOFFUPPER OFFSYMLIST OFFSYMXREF 


OPTIONS 
LIMCOL = 0 , LIMROW = 0 , SOLPRINT = OFF , DECIMALS = 2 
RESLIM = 100 , ITERLIM = 10000 , OPTCR = 0.00 , SEED = 3141; 
k--+----------- DEFINITIONS AND DATA---------------~---------------- 


SCALAR T™ total missiles available /4/ 
AM missiles airborne at one time /2/; 


SETS 
t launch time interval /1*115/ 
tl subset of t /1*25/ 
t2 subset of t /26*50/ 
t3 subset of t /51*75/ 
t4 subset of t /76*100/ 
t5 subset of t /101*115/; 
PARAMETERS 
PK1(t1) negative ln of 1-Pk 
/1 0.2881 
2 0.3173 
3 0.3445 
4 0.3697 
5 0.3927 
6 0.4137 
7 0.4326 
8 0.4495 
9 0.4646 
10 0.4779 
11 0.4897 
12 0.5001 
13 0.5093 
14 0.5173 
15 0.5245 
16 0.5310 
17 0.5368 
18 0.5422 
19 0.5472 
20 0.5520 
21 0.5566 
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22 0.5612 


23 0.5657 
24 0.5702 
25 0.5748/ 
PK2(tZ) negative In of 1-Pk 
/26 0.5794 
27 0.5842 
28 0.5890 
29 0.5939 
30 0.5989 
31 0.6040 
32 0.6091 
33 0.6142 
34 0.6194 
35 0.6245 
36 0.6296 
37 0.6347 
38 0.6397 
39 0.6446 
40 0.6494 
41 0.6541 
42 0.6588 
43 0.6633 
44 0.6677 
45 0.6720 
46 0.6763 
47 0.6804 
48 0.6846 
49 0.6887 
50 0.6928/ 
PK3(t3) negative ln of 1-Pk 
/51 0.6969 
52 0.7011 
53 0.7054 
54 0.7097 
55 0.7142 
56 0.7188 
57 0.7236 
58 0.7285 
59 0.7336 
60 0.7389 
61 0.7445 
62 0.7502 
63 0.7560 
64 0.7621 
65 0.7683 
66 0.7747 
67 0.7811 
68 0.7876 
69 0.7942 


710 0.8007 


0.8072 
0.8136 
0.8199 
0.8259 
0.8317/ 


PK4(t4) negative In of 1-Pk 
0.8373 
0.8425 
0.8474 
0.8519 
0.8560 
0.8597 
0.8631 
0.8661 
0.8688 
0.8711 
0.8733 
0.8753 
0.8771 
0.8789 
0.8808 
0.8828 
0.8850 
0.8874 
0.8901 
0.8930 
0.8961 
0.8993 
0.9023 
0.9050 
0.9070/ 


PK5(t5) negative ln of 1-Pk 

/101 0.9076 

102 0.9063 
0.9021 
0.8941 
0.8812 
0.8618 
0.8348 
0.7985 
0.7518 
0.6935 
0.6229 
0.5396 
0.44393 
0.3364 
0.2181/; 


TABLE 


Cl(t,tl) missile in air confounding matrix 
1234567891011 12 13 1415 1617 18 19 20 21 22 23 24 25 
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4-7 ----------------------MODEL-------------~---- +--+ --------------- 


indicates if firing in time period 


X(t) 
X1(t1) 
X2(t2) 
X3(t3) 
X4(t4) 
X5(t5); 


BINARY VARIABLES 


total pk; 


VARIABLE 
TPK 
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oh 


EQUATIONS 


* MAXIMIZE 
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OBJ.. TPK =E= SUM(t1, X1(t1)*PKI(t1))+SUM(t2, X2(t2)*PK2(t2)) 
+SUM(t3, X3(t3)*PK3(t3))+SUM(t4, X4(t4)*PK4(t4)) 
+SUM(t5, X5(t5)*PK5(t5)); 


* subject to 

NMIS.. SUM(tl, X1(t1))+SUM(t2, X2(t2))+SUM(t3, X3(t3)) 
+SUM(t4, X4(t4))+SUM(t5, X5(t5)) =L= 1™; 

AMIS(t).. SUM(tl, Cl(t,tl)*X1(t1))+SUM(t2, C2(t,t2)*X2(+2)) 
+SUM(t3, C3(t,t3)*X3(t3))+SUM(t4 ,C4(t,t4)*X4(t4) ) 
+SUM(t5 ,C5(t,t5)*x5(t5)) =L= AM; 

ROF(t)..  SUM(tl, R1(t,tl)*X1(t1))+SUM(t2, R2(t,t2)*X2(t2)) 
+SUM(t3, R3(t,t3)*X3(t3))+SUM(t4, R4(t,t4)*X4(t4)) 
+SUM(t5, R5(t,t5)*X5(t5)) =L= 1; 


MODEL MIX /ALL/ 
SOLVE MIX USING MIP MAXIMIZING TPK,; 


* print the optimal objective value and solution. 


DISPLAY TPK.L, X1.L, X2.L, X3.L, X4.L, X5.L; 
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